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Abstract. The colours of the galaxy class known as Extremely Red Objects (EROs; R - K > 5) are considered to be consistent 
with two distinct galaxy populations at high redshift: evolved ellipticals or young dusty starbursts. In this paper the properties 
of EROs, spanned by the five photometric bands RIJHK, are investigated as to the possibility to distinguish between these two 
galaxy populations using only broad band photometry. The broad band colours of elliptical and starburst galaxies at redshifts 
up to 5 are computed from synthetic spectra obtained using the spectral evolution synthesis programme PEGASE.2. Two initial 
mass functions and a range of metallicities and extinctions are used. In order to be extremely red in the redshift range considered, 
the evolution of the R - K colour sets the requirement that ellipticals have to be less than 7-8 Gyr old, and that the starbursts 
must have colour excesses of E(B - V) > 1, as derived from the nebular emission lines. In investigating the overlap in the 
different colour-colour planes as a function of redshift, it is found that the planes formed from permutations of the same three 
filters exhibit very similar overlap characteristics. In colour-colour planes formed within such triplets one of the filters will serve 
as a "pivot" band against which the two other bands are compared. The configuration where this pivot band lies between the 
other two bands has the best performance as a discriminator among the three possible configurations. A consistent behaviour 
cannot be found among the configurations formed by permuting four filters. The minimal filter configuration R - H vs. H - K 
is found to be the very best discriminator, working as such up to redshift 2.9. 

Key words, galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: fundamental parameters - galaxies: high- 
redshift - galaxies: starburst 



1. Introduction 

The optically faint high-redshift galaxy class known as 
Extremely Red O bjects, or EROs, has since its discovery 
(Elst on et al.l ll988) been extensively studied. The EROs are 
considered to be one of the more enigmatic galaxy classes 
since their faintness makes it very difficult to classify them, 
either spectroscopically or morphologically. Although their ex- 
tremely red colours, R - K > 5, narrow down the possibilities 
to essentially ellipticals and dusty starburst galaxies at high 
redshifts, the colour degeneracy between these two very dif- 
ferent galaxy populations prevents the determination of their 
morphologic types. Understanding the nature of the EROs is 
of great importance since they can provide a direct observa- 
tional test of theories of galaxy formati on (e.g. ICimattill200ll): 
the hierarchical cluster ing picture (e.g. Iwhite & Fren kl ll99ll 
iKauffma nn et alj| l993l) versus the mono lifhical collapse pic- 
ture (e.g. lEggen et all 196ilLarsorl 1975b . 

Due to their faintness, only a handful EROs have re- 
liable redshift estimates in combination with a determina- 



tion of the morphologic type. These individual ca ses d ivide 
evenly between starbursts (e.g. ICim atti etal' Ill998t IPey et al. 
ll999HAfonso et al.l200lllsmith et all200lHFraver et al.l2003 ) 



and ellipticals (e.e 
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Send offprint requests to: Stefan Bergstrom, 
e-mail: stefan@oso . Chalmers . se 



conclusion to be drawn as to what galaxy type the majority of 
the EROs belongs to. The individually studied elliptical EROs 
are located in a narrow redshift range, 1.2 < z < 18, whereas 
the individually classified starburst systems span a much larger 
redshift range. 

Surveys have shown that the EROs comprise a heteroge- 
neous galaxy populatio n (Thomps on et al.l 1999- ICimatti etal J 
1999; Moriondo et al. 2000J, although the elliptical case 
has been the stronger one supported by in pa r ticular 
the detected clustering of t he EROs JDaddi et alJ [2000; 
Moustaka s & Somer villel2002l). However, the sample of EROs 
compiled by iDaddi et al J 72002) contains similar fractions 
of star-forming and evolved galaxies, with only the evolved 
galaxies showing strong clustering. In later surveys, the het- 
erogeneity of the ERO population has been confirmed (e.g . 
Man nucci et alJ 120021 ISmail et alJ 12002 1 : ISmifh et alJ |2Q02a; 
Yan & Thompson 2003; Moustakas et al. 2003) and the picture 
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is further complicated by the finding bv lSmith etal1ll2002bl) of 
an early-type disc ERO at z = 1.6. 

Taking the properti es of the ERO sam ple (complete to 
K s < 19.2) obtained by ICimatti et all d2002l) as "typical", the 
ERO population consists of similar fractions of "old" and "star- 
forming" galaxies. The average "old" ERO is concluded to be 
at least 3 Gyr old, being formed beyond z = 2.4, if the metal - 
licity is solar. The average "star-forming" ERO agrees with a 
starburst spectrum having a ICalzettietal.1 d2000h E(B - V) star 
of ~ 0.8, i.e. an extinction for the nebular emission lines of 
E(B - V) = 1.8. If the contribution from an assumed under- 
lying old stellar population is accounted for, the average red- 
dening can be decreased to E(B - V) sta r ~ 0.7. Using other 
synthetic spectra, the extinction is found to be in the range 
0.6 < E(B - V) star < 1.1. 

Since EROs are by their definition very weak at optical 
wavelengths, spectroscopy is not a viable method for classify- 
ing large numbers of EROs. Instead, attem pts have been made 
to use broad band colours. IPozze tti & Ma nnuccil 112000. here- 
after PM00) found that EROs in the redshift ranges 1 < z < 2 
and 2 < z < 2.5 fall in different regions in both / - K vs. J - K 
and R-K\s. J-K plots for sources in the lower redshift range, 
and in H - K vs. J - K plots for sources in the higher redshift 
range. This result is robust and valid for a variety of metal- 
licities, dust contents, star formation histories and redshifts of 
for mation. 

Pieri ni et al .] d2003l) recently extended the study of the 
colour-colour method for classifying EROs. The study is simi- 
lar to that of PM00 but includes dusty post-starburst galaxies as 
well as the nominal dust-free ev olved stellar populations (old 
ellipticals) and dusty starbursts. Pier ini et"al] also studied the 
impact of various dust distributions within the model galaxies. 
While the dust distribution can change the attenuation proper- 
ties, it did not have a large impact on the overall results, which 
were similar to those of PM00. The presence of dusty post- 
starburst systems tends to make the separation of old, dust-free 
systems and dusty starburst systems less distinct. 

In this work, the possibility to use broad band colours to 
distinguish between the elliptical and dusty starburst galaxy 
populations among the EROs at redshifts below five is fur- 
ther pursued, and all possible colour-colour configurations 
of the five broad band filters RIJHK are investigated. The 
study is based on models of the elliptical and dusty starburst 
populations, where synthetic galaxy spectra are obtained us- 
ing the spectral evolution synthesis p rogramme PEGASE.2 
(Fio c & Rocca-Volmerangell997lfl999i) . In particular, colour- 
colour combinations are searched for where the two galaxy 
populations do not overlap over a large redshift range. The in- 
tention is to study to what extent photometric classifications 
of EROs can be used. The adopted cosmology is Hq = 70 
kms 1 Mpc~', Q. m = 0.3 and Qa = 0.7, thus yielding an age of 
the Universe of 13.5 Gyr. 

The paper is organised as follows. The properties of the dif- 
ferent galaxy models are presented in Sect. |2 In Sect. [3] the 
methods used are accounted for, which include determining 
which models produce EROs considering redshift and galaxy 
age, and investigating the overlap between the colours of the 
galaxy populations in different colour-colour diagrams. The re- 



sults are presented and discussed in Sect. @] and Sect. |5] con- 
cludes. 

2. The galaxy models 

In general there are two distinct galaxy populations that the 
EROs are considered to belong to: ellipticals and dusty star- 
bursts. The elliptical galaxies consist of an old stellar popula- 
tion and have none or small amounts of gas and dust. Hence 
the colours of these galaxies are determined by the old stellar 
population which produce a large 4000 A break in the spectral 
energy distribution (SED). The colours of ellipticals can there- 
fore be very red, especially if the galaxies are located at high 
redshift. The other galaxy population producing EROs is the 
dusty starbursts. These galaxies contain much gas and dust and 
have high star formation rates (SFRs). The light of the early 
stellar population is scattered and absorbed by the dust, which 
heats up and subsequently re-radiates the energy in the far- 
infrared. Since this extinction is greater at shorter wavelengths, 
the initially blue SED is reddened. Furthermore, a redshifted 
SED will be redder than the restframe counterpart. 

To model these galaxy population s, the spectral synt hesis 
evolu tion programme PEGASE.2 1 dFioc & Rocca-Vo lmerange 
119971 [ 1999) is used. The programme computes models of 
galactic spectra in the range 91 A to 160 fim assuming dif- 
ferent star formation scenarios, metallicities etc. PEGASE.2 
uses evolutionary tracks of mai nly the Padova group (see 
iFioc & Rocca-Volmerangel Q999 for details), and synthetic 
spectra are used when observational libraries are not avail- 
able. For the nebular emission, only case B recombination is 
taken into account when computing the hydrogen line fluxes. 
For other lines, typical observed HB ratios are used. The pro- 
gramme allows for a consistent metallicity evolution using the 
IWooslev & Weaver! d 1995b models of SNII ejecta. 

The elliptical galaxies are modelled as galaxies that after 
an instantaneous burst of star formation evolve passively, also 
known as a simple stellar population (SSP) mod el. (This is the 
(inst antaneous-)burst population of Bruz ual & Chariot) ll 99 li 
Il993h After the burst of star formation the galaxies are as- 
sumed to contain no interstellar gas or dust, which is equivalent 
to assuming the on-set of a strong galactic wind. The galaxies 
are considered as ellipticals 100 Myr after the starburst, and are 
then let to evolve for a Hubble time. The time grid consists of 
128 epochs with a spacing of about 5% of the galaxy age. 

The starburst galaxies are modelled as galaxies with a con- 
stant star formation during 100 Myr. The time step is 1 Myr. 
The extinction in the starbursts is modelled as a uniform screen 
of dust and gas and follows the ICalzetti et aid lEoOO. hereafter 
COO) obscuration law derived for starbursts. According to the 
work by COO, the stellar continuum and the nebular emission 
lines suffer from different extinctions, due to the different stel- 
lar environments. The extinction of the stellar continuum is 
proportional to the extinction of the nebular emission lines as 
given by the empirically derived relation 



E(B-V) sms =0ME(B - V) gas . 

1 PEGASE.2 is made publicly available at the web site 
http: //www. iap. fr/pegase/. 
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Throughout this paper the extinction will be parameterised us- 
ing the colour excess of the nebular emission lines, and the 
subscript gas is dropped (e.g. COO). 

The COO extinction law expresses the selective absorption, 
not total. The shape of the extinction law depends on the prop- 
erties of the dust grains causing the extinction, whereas the 
amount of extinction can vary significantly. Here, the extinc- 
tion is treated as a screen through which the radiation from the 
galaxies passes, and the amount of extinction can be varied by 
increasing the thickness of the screen. The colour excess re- 
quired to produce an ERO of a particular galaxy at a given red- 
shift is therefore only indicative of the amount of extinction 
needed. The total column density of dust, for example, will 
be dependent on the geometry of the stars and the obscuring 
medium. 

The choice of SFR in the starbursts requires some consid- 
eration as it influences the SED in two ways, (i) The SFR pro- 
vides a scaling factor of the flux of the spectrum. This scaling 
property does not have any impact on the colours as these are 
derived from the shape of the spectrum, not its absolute value. 
This is true also in the presence of nebular emission lines, as the 
continuum flux and emission line flux (through the flux of ion- 
ising photons) both are proportional to the number of stars pro- 
duced. This has the result that two galaxies of the same metal- 
licity, IMF and temporal shape of the star formation history will 
have the same colours regardless of the magnitude of the SFR. 
(ii) The SFR will influence the rate at which the galaxies are 
enriched: the higher the SFR, the faster the enrichment. This 
will cause a major change in the colours. Hence, of these two 
effects it is only the latter, where the metallicity is increased, 
that has a significant effect on the colours. 

Since a consistent evolution of the metallicity is used in the 
modelling of the galaxy spectra, the metallicities will change 
in a galaxy experiencing an extended period of star forma- 
tion (i.e. the starburst model). Following the line of reasoning 
sketched out in the previous paragraph, the SFRs cannot be too 
high if the range of metallicities is to be kept in the modelling. 
Hence, the value of the SFR is chosen to achieve a well con- 
trolled range of metallicities. The metallicity grid consists of 
five metallicities, 0.02, 0.2, 0.4, 1 and 2.5 Z Q , which in the pas- 
sive evolution of the ellipticals do not change, but in the case 
of the starbursts, these are initial values. 

The SFR chosen to ensure such a slow enrichment is 
5 x 10 _1 'Mgaiaxy yr _1 . A galaxy with such an SFR would not be 
regarded as a starburst galaxy since only 0.5% of the gas mass 
have been used to form stars after 100 Myr. However, since the 
magnitude of the SFR does not influence the colours of galax- 
ies of the same metallicity and stellar population, as explained 
above, this low-SFR model is equivalent to models having a 
more "realistic" SFR. Considering the metallicity in the ISM, 
an initial solar metallicity has been increased to 1 .0050 Z , i.e. 
by 0.5%, after 100 Myr for a starburst of Salpeter IMF. The 
enrichment is stronger the lower the initial metallicity, and the 
highest increase of metallicity occurs for the starbursts of ini- 
tial metallicity 0.02 Z G and Salpeter IMF; after 100 Myr, the 
metallicity have increased by ~25%. 

Th e initial mass fu nctions (IMFs) used in the mo delling 
are the ISalpeteJ \ 19551) and the "present-day" iKrouoal J200ll 
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Fig. 1. Spectral energy distributions of galaxies at z — 2 having 
Salpeter IMF and solar metallicity, normalised in the K band. 
The solid line is the SED of a 1 Gyr old elliptical galaxy; the 
dashed line is the SED of a 100 Myr old starburst galaxy having 
E(B — V) = 2. Shown are also the wavelength regions covered 
by the R, H and K filters. 



Table 1. PEGASE.2 parameters common to all the models 



Property 


Value 


Supernova ejecta 


Wooslev & Weaver (1995) 




model B 


Infall 


No 


Stellar winds 


Yes 


Fraction of close 




binary systems 


0.05 


Metallicity evolution 


Consistent 


Formation of substellar 




objects 


No 


Nebular emission 


Yes 


Extinction 


No 



equation (6)) IMFs in the mass range 0.09 < M/M < 120. 
The Kr oupa I MF is mainly based on the power-law compila- 
tion by Scalo ( 1998) but the effects of unresolved binaries are 
also taken into account. This IMF differs from the Salpeter IMF 
below 1 M Q and is included to increase the generality of the re- 
sults. The adopted mass limits are consistent with the stellar 
tracks used by PEGASE.2. The parameters common to all the 
models are shown in Tabled 

Model spectra of a 1 Gyr old elliptical and a 100 Myr old 
starburst galaxy of E(B-V) = 2 are shown in Fig.^ The galax- 
ies have Salpeter IMF and solar metallicity. Inlaid are the wave- 
length regions of the R, H and K bands. In Fig. [2] spectra of 
galaxies of ellipticals 200, 1000 and 5000 Myr old (left panel) 
and of starbursts 0, 10 and 100 Myr old having E(B -V) = 2 
(right panel) are shown. 
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Fig. 2. Spectral energy distributions of galaxies having Salpeter IMF and solar metallicity, normalised in the K band. Left panel: 
a 200 Myr, 1 Gyr and 5 Gyr old elliptical. Right panel: a new-born, 10 and 100 Myr old starburst galaxy having E(B — V) = 2. 
Both galaxy types become redder with age. 



3. The method 

Two questions are addressed in the present work: (i) When can 
ellipticals and starbursts be EROs? (ii) Is it possible to distin- 
guish between the extremely red ellipticals and dusty starbursts 
merely from their broad band colo urs? A first attempt to an swer 
the second question was made by Scodeggi o & SilvalfeOOOl) . In 
a similar study, PM00 were able to separate old elliptical galax- 
ies from dusty starbursts using / - K vs. J - K and R - K vs. 
J - K plots in the redshift range 1 < z < 2, and also H - K vs. 
J - K in the range 2 < z < 2.5. 

In the work presented here, the ability to distinguish be- 
tween the considered galaxy types are investigated for all the 
45 possible combinations of the RIJHK filters and for all plau- 
sible redshift ranges. This section describes the methods used 
to investigate the colour evolution and the colour plane proper- 
ties. 

The fluxes at different redshifts are computed from the 
spectra obtained from PEGASE.2 using the standard relation 



F(Ao) 



LQ*i)(l+z) 
AnD\ ' 



(2) 



where Aq = (1 + z)M and Dl is the luminosity distance. The 
magnitudes are then computed in the Vega system using the 
well-documen ted filter functions of B essell & Brettl |l988 ) and 
iBesselll |l990) founded on the Cousins system. The filter re- 
sponse functions are extended to 1500 data points using cubic 
splines. This filter system is very similar to the filters used in 
the EIS-Deep Survey, the SUSI RI and SOFI JHK filters at the 
NTT 2 , so any conclusions made in the following will not be 
altered significantly if the latter system is used. 



Table 2. The properties of the normalised (unity peak transmis- 
sion) Bessell filters. The width is defined as the area of the filter 
response function. 



Band 


^mean (A) 


Width (A) 


F (Wnr 2 A" 1 ) 


R 


6586 


1581 


2.216 x 10~ 12 


I 


8060 


1495 


1.154 x 10 12 


.1 


12368 


2023 


3.256 x 10- 13 


H 


16466 


2855 


1.181 x 10~ 13 


K 


22119 


3660 


4.110 x 10~ 14 



2 Filter functions achieved from the Filter Catalogue at the ESO La 
Silla home page, http: //www. Is . eso . org. 



3.1. Colour evolution 

Two approaches are used to find galaxies that might produce 
EROs under certain circumstances. The first approach is to find 
the lowest redshift required to make a galaxy of a certain age 
extremely red. The redshift where R - K — 5 will hereafter be 
called zero- This procedure imposes a cosmology consistency 
criterion: the redshift of formation, zi, of a galaxy observed (as 
an ERO for instance) at a certain redshift has to be consistent 
with the assumed cosmology. As will be seen, this cosmologi- 
cal constraint rules out that ellipticals capable of being classi- 
fied as EROs are older than about 7-8 Gyr. 

The second approach is to study how the R - K colour 
evolves with redshift. This can be done both by assuming a 
redshift of formation, zu and let a galaxy evolve in time, or by 
redshifting a galaxy of a fixed age. 



3.2. Overlap in the colour plane 

The main aim of this work is to investigate whether broad band 
colour-colour plots can be used to discriminate between the 
two main classes of galaxies making up the ERO population: 
evolved ellipticals and dusty starbursts. The colours of ellip- 
tical and starburst galaxies at different redshifts may or may 
not overlap, and there could exist redshift ranges in which the 
colour-colour regions of the two galaxy populations can be dis- 
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tinctly defined. To investigate this, the first step is to build up a 
grid of colours of ERO-classified galaxies. For a range of red- 
shifts, the colours of ellipticals and starbursts are computed. 
After removing those that are neither consistent with the cos- 
mology nor have R-K > 5, two sets of colour-colour points are 
obtained, associated with the EROs in each galaxy population. 

The overlap that occurs in a colour-colour diagram between 
the colours of the two galaxy populations at different redshifts 
is then computed. The confusion, or overlap, between these two 
point sets are defined as follows. If A designates the region cov- 
ered by the one set, and B the region of the other set, then the 
intersection of these two sets, A n B, is the region where the 
two sets overlap. The measure of overlap is now taken as the 
area of the intersection normalised by the area of the smallest 
region, i.e. 

An B iAnB AnB\ 

overlap = = max { , } . (3) 

F min{A,B) I A B I 

In this way, an overlap measure of 1 is obtained if one of the 
sets is completely embedded in the other. 

To define the regions covered by the point sets, every 
colour-colour point is considered as defining the centre of a 
circle having a fixed (small) radius. If the point density is high 
enough, this will create a continuous region. Two points are 
now regarded as overlapping if the distance between them is 
less than twice the radii of the circles. This requires that the 
point density is homogeneous, since it cannot be accepted that 
a higher point density part of the colour-colour region will have 
a higher weighting in the overlap computation. Hence, the two 
point sets are discretised into an evenly spaced grid, with a 
spacing of one circle radius. The colour-colour points will be 
offset from their true loci by at most one half of a circle radius. 
In the overlap computation, the circle radius is set to 0.05 mag- 
nitudes, which will make two points more than 0.1 magnitudes 
apart seen as separate. The gridding will shift the colour-colour 
points by at most 0.025 magnitudes, which is well below the 
typical photometric errors. 

This procedure results in maps of overlap measures where 
one of the axes corresponds to the redshifts of the ellipticals, 
and the other to the redshifts of the starbursts (see Fig. [6). 
From these maps, the redshift regions where the starburst EROs 
and the elliptical EROs can be separated from each other in a 
colour-colour diagram can be clearly seen. What is left is to 
explore these redshift regions to see what colours the different 
galaxies have. 

4. Results and discussion 

4.1. Colour evolution 
4.1 .1 . The ellipticals 

All ellipticals have R - K > 5 at some redshift and the general 
trend of the R-K colour is to increase with redshift. The expla- 
nation for this can be seen in the left panel of Fig. [3] when the 
SED of an elliptical is redshifted, it will be shifted to the right in 
the diagram, and except for the youngest ellipticals, the R filter 
will soon measure in the less luminous region blue-ward of the 
4000 A break, while the K filter measures brighter and brighter 




i 1 1 1 1 1 1 1 1 

2 4 6 8 

Galaxy age (Gyr) 

Fig. 3. zero vs. galaxy age for the ellipticals of Salpeter IMF, 
where zero is defined as the lowest redshift where an ellipti- 
cal can be an ERO. The metallicities (0.02 - 2.5 Z Q ) of the 
models increase downwards/leftwards and the lines thicken ac- 
cordingly. The dashed line traces the cosmological redshift-age 
relation. 



parts of the SED. A higher metallicity reddens the SED due to 
line-blanketing and backwarming in the stellar atmospheres. A 
redder galaxy will require a lower redshift to make the galaxy 
extremely red, i.e. a lower zero (as defined in Sect. I3.U . As a 
result of this, an elliptical of higher metallicity can be encoun- 
tered as an ERO over a larger redshift range and hence it can 
also be older compared to an ERO elliptical of lower metallic- 
ity. 

In Fig. [3] the Zero s f° r the ellipticals are plotted as a 
function of galaxy age. The dashed line shows the redshift- 
age correspondence for the chosen cosmology: the area left- 
wards/downwards of this line contains the redshift-age set con- 
sistent with the cosmology. The five solid lines show the zeros 
for ellipticals of different metallicities (0.02 - 2.5 Z Q ), and 
the thickness of the lines increases with increasing metallicity. 
Since zero is the lowest redshift needed to produce an ERO, the 
redshift regime between the "cosmology line" and each zero 
line shows the redshifts where the EROs exist. It is clearly seen 
that the cosmology requires the ellipticals to be younger than 
7.5 Gyr. Furthermore, the lowest zero is abo ut 0.7, a r esult 
which is consistent with the observations of Moustakas et all 
(2003), who find very few early-type EROs at redshifts below 
-0.9. 

The elliptical R-K colour for different assumed redshifts 
of formation is shown in Fig. 0] The lowest metallicity, Z = 
0.02 Z , only produces EROs in a small redshift interval around 
z = 2 for the highest Zf. For decreasing Zf, the limiting ERO- 
producing metallicity increases, and when Zf = 1, not even the 
highest metallicity will yield red enough galaxies. 
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012345 012345 
Redshift Redshift 

Fig. 4. Colour-redshift diagrams for the ellipticals of Salpeter IMF assuming formation redshifts from 2 to 5 using the full range 
of metallicities (0.02 - 2.5 Z ). The thickness of the lines increases according to the metallicity. 



4.1.2. The starbursts 

For starburst galaxies, the evolution of the R - K colour with 
redshift puts constraints on the amount of reddening needed 
to produce EROs. This limiting amount of extinction will de- 
crease as the metallicity increases, since a more metal-rich 
galaxy will be redder. The colour-redshift diagram in Fig. |5] 
shows the regions of the R-K colour of the starbursts for colour 
excesses of 1 .0 and 2.0. The full ranges of galaxy age (0 - 100 
Myr) and metallicity (0.02 - 2.5 Z Q ) are used. A colour ex- 
cess of slightly more than 1 is required to produce EROs for 
Z < 5. Excepting the very lowest redshifts, the starbursts can be 
observed as EROs at all z < 5 when E(B - V) = 2.0. 

4.2. The colour plane properties 

The following investigation of the colour-colour space, both as 
the overlap maps and the colour-colour plots are concerned, 
does not depend on any assumed formation redshifts. The 
galaxies are assumed to form continuously, and only the cos- 
mology sets the age limits for the galaxies. At a certain redshift 
the considered galaxy models are those that strictly fit into the 
current lifetime of the Universe. No assumption is made of a 



1 1 1 1 1 
^^E(B-V) = 1.0 


i i i i i i i 

















I I I I I I I I I I I I I I I I I I I I I I I I I I 

1 2 3 4 5 



Redshift 

Fig. 5. Colour-redshift diagram for the starbursts of Salpeter 
IMF. The vertically hatched region corresponds to E(B - V) = 
2.0, and the horizontally hatched region to E(B — V) = 1 .0. The 
full age and metallicity ranges are used; the more metal-rich 
and the older the galaxy, the redder. 
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redshift where the galaxy formation is started in the Universe 
as a whole. 

4.2.1 . The overlap in the colour-colour diagrams 

Distinguishing between the elliptical and starburst models in 
colour-colour diagrams is possible only if the colours of the 
galaxies follow evolutionary tracks that do not overlap. As a 
first test, colour-colour plots of all the 45 colour-colour config- 
urations that can be formed from the filters RIJHK are made 
for the whole redshift grid (redshifts below 5 in steps of 0.1). 
Both Salpeter and Kroupa IMFs are used, and the colour ex- 
cess of the starbursts ranges up to 2.0 in steps of 0. 1 . It is found 
that when considering this broad redshift range, all filter con- 
figurations produce colour-colour plots where the region of the 
starburst colours lie almost completely inside the colour region 
of the ellipticals, hence they cannot be separated in this type of 
plots. 

The next step is to see in which redshift ranges the different 
configurations produce colours where the galaxy populations 
are distinguishable, and the overlap is computed using the pro- 
cedure described in Sect. 13.21 The model grid is again set up 
using Salpeter and Kroupa IMFs and colour excesses of the 
starbursts up to 2 in steps of 0.1. The redshift grid spacing is 
0.1. Overlap maps resulting from this procedure are shown in 
Fig.0 

A proper interpretation of these maps requires some re- 
marks, (i) The overlap has to be small over large continuous 
regions in the redshift space, starting from z — and up to 
high redshifts. A colour configuration is of course only use- 
ful as a discriminator at redshifts below the redshift where 
the colours starts to overlap, (ii) Since the fraction of EROs 
in a given galaxy population increases with redshift, the re- 
gion occupied by EROs in a colour-colour plot will increase 
with redshift as well. At low redshifts, the colour region cov- 
ered by the starbursts can therefore be very small, which will 
make any overlap with the ellipticals to appear more significant 
than what is reasonable. Hence it is not enough to decide on a 
redshift limit for a certain filter configuration from its overlap 
map, but to determine the impact of the overlap at low red- 
shifts, the corresponding colour-colour diagram has to be in- 
vestigated, (iii) The ridges in the overlap maps often represent 
cases where the colours of one galaxy population cross over the 
colour region of the other. This causes the colour-colour plots 
to exhibit very complicated structures which makes it difficult 
to define a separation line. 

The configurations that can be formed from permutations 
of a set of three bands (e.g. the triplet {RHK} from which 
the filter configurations RHHK 3 , RKHK and RHRK may be 
formed) have a "pivot band" against which the two other mag- 
nitudes are compared, e.g. the K band in the RKHK configu- 
ration. The configurations formed from such a triplet are found 
to have very similar overlap features, although small differ- 
ences can be seen. For all the triplets, the configurations where 
the pivot wavelength band lies between the other two bands, 

3 The notation "XYZW" will be used as an abbreviated form of "X— 
Yvs.Z- vr. 



are always the best ones; for instance RHHK in the triplet 
{RHK}. It should be noted that this result is consistently ob- 
served throughout the course of this work. The worst perfor- 
mance is in most cases exhibited by the configurations us- 
ing the shorter wavelength band as pivot band, although some 
(small) changes can be seen if higher extinctions are included. 

The configurations formed from quadruples (configurations 
made with four different bands, e.g. the quadruple {RJHK} 
from which the configurations RJHK, RHJK and RKJH may 
be formed) tend to behave very differently and cannot be re- 
garded as forming distinct groups having similar overlap fea- 
tures. Since they are based on four regions on the SED, they 
are more sensitive to differences between the SEDs. However, 
in the overlap maps this advantage over the triplets manifests 
itself only in that the average performance of the quadruples is 
better than the triplets, but the best triplet has still a higher red- 
shift limit than the best quadruple configuration. Furthermore, 
the only quadruple that actually performs better than the four 
triplets it can be broken into is the {IJHK} quadruple. This 
combination of filters appears unfortunate; the triplets {IHK}, 
[JHK] and {IJH} are the three worst of the triplets, and the 
{UK} as well as [IJHK] itself both yield complicated colour- 
colour plots. For the other quadruples, there is always a con- 
figuration of three filters that performs at least as well as its 
"parent" quadruple. 

The relative discriminating performance of the three con- 
figurations within each filter set are shown in Table [3] Figure 
|6]shows the filter configuration yielding the best overlap prop- 
erties within a given filter set. The colour coding means the 
brighter the area, the less overlap. The RHHK configuration 
has one of the highest redshift limits, and can discriminate be- 
tween the galaxy populations up to redshift 2.9. The small over- 
lap between ellipticals and starbursts at z ~ 1 for this configu- 
ration does not pose a problem. This will be discussed further 
in the next section. 

4.2.2. Non-overlapping colour regions 

Having found the redshift regimes where the best colour-colour 
plots exhibit distinct elliptical and starburst colour regions as 
described in Sect. l4.2.T1 an attempt is now made to define lines 
that separate the two regions. A filter configuration can be con- 
sidered as having an improved performance as a discriminator 
if the separation line is parallel with the reddening vector, since 
an increase in extinction will not cause the starbursts to cross 
such a separation line. However, as a higher extinction also will 
enable lower-redshift starbursts to become EROs, it is possible 
that these "new" candidates will intrude on the regions of the 
ellipticals. 

A separation line between the colour regions of ellipticals 
and starbursts are defined for the configuration having the best 
performance within each filter set, see Table [3] For a given 
colour-colour configuration, it is in many cases possible to de- 
fine a whole set of lines separating the two regions; in these 
cases a representative line is presented. The separation line is 
valid for z < Zmax, where z max is the redshift limit beyond which 
the colours will be confused. A separation line is not defined for 
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Fig. 6. The colour overlap maps for the filter configurations yielding the best discriminator performance within each filter set. The 
brighter the area, the less overlap; white corresponds to a separation of the elliptical and starburst EROs by at least 0.1 magnitudes 
in the corresponding colour-colour plot. The galaxy models have Salpeter and Kroupa IMFs, metallicities from 0.02 to 2.5 Z , 
and for the starbursts, colour excesses up to 2.0. Only models yielding EROs consistent with the cosmology are considered in 
computing the colour overlap. 
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Table 3. This table shows (i) The relative performance of the filter configurations within each filter set. This is based on the 
ability to discriminate between the colours of the modelled ellipticals and starbursts for E(B — V) < 2. (ii) The separation line 
corresponding to the configuration that has the best discriminating performance within each filter set. (iii) The redshift limit 
beyond which the colour regions become entangled, z max . The separation line is valid up to this redshift limit. Separation lines 
are not defined for the filter configurations having a redshift limit less than about one. (iv) The faintest ERO K magnitude of L* 
ellipticals and starbursts at z max . 



Filter set 


Configuration performance 


Separating 


line: colours of 


vmax 




L*) at Zmzx 




Best 


— * 


Worst 


ERO ellipticals at z < z m ax 




Ellipticals 


Starbursts 


{RHK) 


RHHK 


RKHK 


RHRK 


(R-H)> 


3.69 (H-K) + 0.60 


2.9 


20.7 


26.1 


{RJK} 


RJJK 


RKJK 


RJRK 


(R-J) > 


1.99 (J - K) - 1.18 a 


2.2 


20.7 


25.4 


{UK} 


IJJK 


IKJK 


IJIK 


(I -J) > 


1.21 (J - K) -0.64 ab 


2.2 


20.7 


25.4 


{RJH\ 


RJJH 


RHJH 


RJRH 


(R-J) > 


1.77 (J-H) + 1.10 


2.1 


20.7 


25.3 


{RIJ} 


RIIJ 


RJIJ 


RIRJ 


(R-I) > 


-0.38 (I -J) +2.05 


1.6 


20.7 


24.8 


IRIK] 


RIIK 


RIRK 


RKIK 


(R-I) > 


0.22 (I-K) +0.39 


1.4 


20.6 


24.6 


{RIH} 


RI1H 


RHIH 


RIRH 


(R-I) > 


0.25 (I-H) +0.56 


1.2 


20.5 


24.4 


{IHK} 


IHHK 


IKHK 


IHIK 












{JHK} 


JHHK 


JKHK 


JHJK 












[IJH] 


UJH 


IHJH 


IJIH 












{RJHK} 


RJHK 


RHJK 


RKJH 


(R-J) > 


4.59 (H-K)- 1.44 a 


2.6 


20.7 


25.8 


{1JHK) 


IHJK 


IJHK 


IKJH 


(I-H) > 


1.78 (J - K) -0.72 ab 


2.4 


20.7 


25.6 


{RIHK} 


RIHK 


RH1K 


RKIH 


(R-I) > 


1.80 (H - K) - 0.59" 


2.0 


20.7 


25.2 


{RIJK] 


RIJK 


RJIK 


RKIJ 


(R-I) > 


0.78 (J - K) -0.35" 


1.8 


20.7 


25.0 


{RIJH} 


RJIH 


RIJH 


RH1J 


(R-J) > 


0.65 (I-H) + 1.20 


1.8 


20.7 


25.0 



11 Line parallel with reddening vector. 

b Line runs partly through the colour-colour regions due to complex region shapes. 



Table 4. The consistency of the classifications made using different filter configurations are checked against classifications of 
samples of ER Os in the liter ature. The consistency is grouped accor ding to classification method: (i) the photometric classifi- 
cation s bvlSmail et al. d2002l); (ii) the spectroscopic classificat i ons by [Graham & Devl(ll996T).ISpinrad et al ] Jl997h . lSoifer etalJ 
Jl999l) . lAfonso et alJ feOOlUSmith et al. (l200ll).lFraver et alllEooi) andlSaraccoetal.ld2003l): and (iiiVthe morphologic clas- 
sifications bv iFranceschini et alJ Jl998t) . IStiavelli et all i 1 9991) . ISmith et alJ J2002al)~ d Ivan Dok kum & Stanford! J2003I) . The 
consistency is given both as a colour consistency, i.e. the consistency with the separation lines given in Table and a region 
consistency, i.e. if the sources fall in the proper colour-colour regions. The consistency is given as the number of sources on the 
correct side of the separation line (or in the correct colour region) over the respective number of sources available. 



Filter set 


Configuration 


Zmax 


Colour (region) consistency for different classification methods 








Smail et al. 2002 sample 


Spectroscopically classified 


Morphologically classified 








Evolved 


Dusty 


Evolved 


Dusty 


Evolved 


Dusty 


{RHK) 


RHHK 


2.9 


6(6) / 8 


11(11) / 15 


2(2) / 2 


1(1)/ 1 


1(1)/ 1 




{RJK} 


RJJK 


2.2 


9(8) / 9 


3(1") / 9 


4(4) / 4 


2(2) / 3 


3(3) / 4 


0(0) / 10 


{UK} 


IJJK 


2.2 


7(4) / 7 


4(3) / 10 


3(2) / 3 


0(0 a ) / 3 


20(11) / 21 


1(1")/ 9 


{RJH} 


RJJH 


2.1 


4(2) / 5 


KD/ 9 


2(2) / 2 


0(0) / 1 


0(0) / 1 




{RIJ} 


RIIJ 


1.6 


5(3) / 7 


5(3) / 9 


3(1)/ 3 


1(0) / 2 


4(2) / 4 


0(0) / 9 


{RIK} 


RIIK 


1.4 


3(2) / 10 


11(10) / 17 


KD/ 3 


l(0 a ) / 2 


4(4) / 8 


9(7 a ) / 20 


{RIH} 


RIIH 


1.2 


2(0) / 8 


9(6) / 15 


1(0) / 1 


0(0) / 1 


0(0) / 1 




{RJHK} 


RJHK 


2.6 


5(5) / 5 


7(3 a ) / 9 


2(2) / 2 


K0 a ) / 1 


KD / 1 




{IJHK} 


IHJK 


2.4 


5(4) / 5 


5(3) / 10 


KD / 1 


1(1)/ 2 


17(11) / 18 




{RIHK} 


RIHK 


2.0 


5(5) / 8 


ll(7 a ) / 15 


KD / 1 


l(0 a ) / 1 


0(0) / 1 




{RIJK} 


RIJK 


1.8 


5(3) / 7 


6(3) / 9 


3(3) / 3 


2(0 a ) / 2 


2(2) / 4 


4(3 a ) / 9 


{RIJH} 


RJIH 


1.8 


4(2) / 5 


2(0) / 9 


KD / 1 


0(0) / 1 


0(0) / 1 





a The region consistency can be improved by using a higher extinction. 



the filter sets that have redshift limits below about one, as non- 
overlap at such low redshifts are due to a lack of EROs among 
the ellipticals and/or the starbursts. 

To check the performance of the different filter configu- 
rations, comparisons are made with existing classifications in 



the literature. This is shown in Table |4] The consistency of 
the classifications made here is checked according to classi- 
fication method: photometric, spectroscopic or morphologic. 
The photometri c cons istency check is made against the data 
in ISmail et alJ 00021 S02 hereafter). S02 have obtained a 
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sample of EROs from deep RIzJHK photometry and have 
also classified the sources by template fitting. The sample is 
thus divided into two classes where the photometry is con- 
sistent with "dusty" or "evolved" galaxi es. The spectrocop - 



Spinrad et al 
Smith et alJ 




SoiferetalJJl999). 


Afonso et al. (2001), 


Fraveretal.1 12003 


) and 


Saracco et alJ 



2003 ) . The morphologically clas si fied EROs are ob - 
tained from iFranceschini et al.1 rtl998h. IStiavelli et alJ |l999), 
Smith et alJ J2002al> and Ivan Dokkum & Stanford! J2003). 
Smi th et alJ ( I2002ah use the classification "Irregular" and 
"Compact", which, for the consistency check, is simply 
interpreted as "starburst" and "elliptical" (or "dusty" and 
"evolved"). This might be too simplistic which would explain 
the apparently large discre pancy for the morphologically clas- 
sified starbursts, since the Smi th et alJ (l2002al) sample com- 
prises the majority of these sources. 

The RHHK configuration is the best one as it concerns red- 
shift, and it also agrees well with existing classifications. Figure 
[TJshows the RHHK colour diagram for starbursts and ellipticals 
up to redshift 2.9. This filter configuration agrees with all the 
spectroscopically and morphologically classified sources - al- 
though the number of sources with H band data is limited - 
and with the classifications made by S02 for 11 of 15 dusty 
and 6 of 8 evolved galaxies that have data in these three bands. 
Furthermore, one of the 1 1 unclassified S02 sources is found to 
have the colours of starbursts and 8 (possibly 10) resemble el- 
lipticals. The sources that end up inside the "wrong" colour re- 
gion can be explained as lower-redshift starbursts having colour 
excesses between 2 and 3 and ellipticals at redshifts higher than 
2.9, respectively. This combination has the property of being a 
minimal choice of filters, since it only requires one extra fil- 
ter except for the ERO criterion filters. Figure^shows how the 
R -K colour can select the extremely red elliptical and starburst 
galaxies, and how they can be distinguished between in the H 
band. 

The RJJK configuration has a smaller redshift coverage, 
but is also a minimal choice of filters. However, this does not 
agree with the existing classifications very well, in particular 
the consistency for the starbursts is bad. In this configuration, a 
majority of the sources have colours that are consistent with the 
colours of the model ellipticals. The mismatch also occurs for 
the {UK} configurations, where additionally the complicated 
structure of the colour-colour regions makes it difficult to de- 
fine a clear separation line. The configurations in the {IHK}, 
{JHK} and [UH] triplets have all very high overlap at low 
redshifts, which can be attributed to the problem lined out in 
remark (ii) in Sect. 14. 2.11 above. However, as these high over- 
lap regions extends over such considerable redshift space, these 
filter configurations must be regarded as having a poor discrim- 
inating ability. 

The best quadruple configuration is RJHK. The separation 
line in this case agrees well with earlier classifications, but a 
majority of the "dusty" S02 sources fall slightly outside the 
starburst region. A heavier extinction could justify this, though. 

A problem with this kind of method to classify EROs is in- 
dicated by the last two columns in Table [3] These show the K 
magnitude of the faintest L* (2 x 10 10 L ) ERO at the limit- 



I CD 




H K 

Fig. 7. R - H vs. H - K regions for the Salpeter and Kroupa 
IMF models for z < 2.9 having extinctions E(B - V) < 2, 
plotted together with classified EROs (references in text and in 
Tabled. Stars denote dusty galaxies, circles evolv ed galaxies, 
and cr osses denote the unclassified detections by ISmail et alJ 
The arrow indicates how the starburst region shifts if 
the colour excess is increased by 0.1. The error bar shows the 
typical photometric error of the colours of the EROs. 



ing redshift for each colour plane. In general, an L* starburst 
ERO is 4 - 6 magnitudes fainter in the K band than an ellip- 
tical at z max . This points to a severe difficulty for studies of 
EROs aimed at disentangling this population: when a study 
is deep enough to detect "ordinary" starburst galaxies, ellipti- 
cals at redshifts far beyond z max will also be detected. Without 
redshift information, these high-redshift ellipticals will be in- 
cluded in the sample, and the result will be colour planes with 
a large amount of confusion. On the other hand, a more shal- 
low survey will tend to detect a majority of ellipticals, which 
can mislead the conclusions. 

In its extension, this problem can affect conclusions drawn 
from studies aimed at determining the fractions of different 
galaxy types among the EROs using any of the standard meth- 
ods. Considering these differences in magnitude of the two 
populations, it seems plausible that the conclusion that ellipti- 
cals dominate the ERO population might just be an effect of this 
bias. However, a thorough study of the luminosity functions of 
the galaxy populations constituting the ERO class is required to 
resolve these implications. The widely different characteristics 
of the two main candidate populations will presumably make 
studies aimed at determining the relative fractions of galaxy 
types among the EROs very sensitive to biases. 

4.2.3. Redder EROs 

The effect of applying higher extinctions is that a greater num- 
ber of lower-redshift starbursts will be included in the ERO 
population, some of which will have colours similar to the el- 
lipticals. As mentioned before, the starbursts which are on the 
correct side of a separation line that is parallel to the redden- 
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J-K 

Fig. 8. R - 7 vs. J-K regions for the Salpeter and Kroupa 
IMF models for z < 2 having extinctions E(B - V) < 3, plot- 
ted together with classified EROs (references in text and in 
Table|4}. Stars denote dusty galaxies, circles evolv ed galaxies, 
and cr osses denote the unclassified detections by ISmail et alJ 
ll2002l) . The arrow indicates how the starburst region shifts if 
the colour excess is increased by 0.1. The error bar shows the 
typical photometric error of the colours of the EROs. 



ing vector will not cross the line if the extinction is increased. 
However, the new ERO candidates among the low-redshift star- 
bursts which previously had R-K < 5 can, as they are reddened 
further, appear "from below" and on the wrong side of the sep- 
aration line, thus ending up in the region of the ellipticals. A 
conclusion from this, is that the redshift plays a major part in 
determining which side of the separation line that the starbursts 
populate. In fact, this appears to be a general effect of the red- 
shift, since if higher-redshift ellipticals are included, these will 
move towards the separation line as the redshift is increased, 
and eventually cross the line. 

A higher extinction does in general affect the triplet con- 
figurations more than the quadruples. For E(B — V) < 3, the 
most stable triplet configuration is RJJK, which is valid for 
Z < 2, although it has a considerable overlap between starbursts 
at z < 0.5 and ellipticals at z ~ 0.8. Accepting this overlap, the 
proper separation line is (R - J) = 1.99(7 - K) — 0.70; see 
Fig. The best quadruple configuration when E(B — V) < 3 
is RJHK, in which the ellipticals have colours consistent with 
(R-J)> 0.90(7 - K) - 0.28 for z < 1.8. Note that these two 
separation lines are parallel with the reddening vector. 

The increased overlap at low redshifts caused by high ex- 
tinction starbursts can be remedied by using a more strict ERO 
criterion, such as R - K > 6, which will disentangle the colours 
by simply de-selecting galaxies, in particular those at z 55 1, cf. 
Figs. |4] and [5] Hence, a stricter criterion will stabilise the dis- 
criminating properties of the colour planes, and also in general 
increase their redshift limits. When the criterion R — K > 6 is 
used, the two configurations that have the best redshift cover- 
age are the RJJK (ellipticals have (R - J) > 1 .99(7 - K) - 0.94 
for z < 2.1) and the IJJK ((/ - J) > 1.21(7 - K) - 0.30 for 



z < 2.2). Also in these cases, the separation lines are parallel 
to the reddening vector. The latter configuration will require an 
ERO selection criterion, though. 

The RHHK turns out not to be very stable for an increase 
of the starburst extinction when using the standard ERO cri- 
terion (R - K > 5), as lower-redshift starbursts of high ex- 
tinction are selected and fall in the colour region of the ellip- 
ticals. Using the stricter ERO criterion R - K > 6, RHHK 
is the best configuration for E(B — V) < 2 and is now valid 
for z < 3.2. High-redshift ERO ellipticals now forces the sep- 
aration line to the right in Fig. [7J these are consistent with 
(R - H) > 4.00(77 - K) - 0.15. Pushing the colour excess to 
higher values shows that RHHK can be used for E(B -V) < 2.4 
when z < 3. For this high extinction, the position of the separa- 
tion line is dictated by the starbursts, and is now moved slightly 
to the left, to (R - H) > 3.69(77 - K) + 0.77. 

Due to a possible shortcoming of existing spectral evolution 
synthesis programmes, the modelled SEDs might in fact be too 
blue. As pointed out bv lZackrisson et al . (2001), the pre-main 
sequence evolution of stars will have a non-negligible effect 
on the SED during the first tens of Myr for a continuous star 
formation, a fact which is usually overlooked in existing pro- 
grammes, such as PEGASE.2, among others. However, as this 
phase reddens the SED, this is not a serious problem for this 
study, because (i) in colour-colour diagrams, such as Figs. [7] 
and the starbursts will move up and to the right, and hence 
away from the loci of the ellipticals, and (ii) since this is a tran- 
sient effect, the younger starbursts are moved closer to the older 
starbursts in Fig. [5] which actually narrows the area covered, 
hence working to decrease the overlap. 

5. Summary and conclusions 

In this paper the colour-colour properties of EROs are inves- 
tigated. Two types of galaxies, starbursts and ellipticals, have 
been modelled using the spectral synthesis evolution program 
PEGASE.2. The ellipticals are modelled as dust-free galaxies, 
passively evolving after an instantaneous burst of star forma- 
tion. The starbursts are modelled as having a constant star for- 
mation rate for 100 Myr, with dust extinction of variou s degrees 
follow ing the Calzetti et al. (2000) obscuration law. ISalpeterl 
( 1955) and Kroupa (2001, equation (6)) IMFs are used and the 
metallicities of the galaxies were evolved consistently with age, 
using initial metallicities from 0.02 up to 2.5 Z . The colour 
properties in the 45 colour planes formed from the RIJHK 
bands are explored up to redshift 5 in a 77o = 70 km s _1 Mpc~', 
Q m = 0.3 and Qa = 0.7 cosmology. 

The following conclusions can be drawn. 

1 . Except for the very lowest-metallicity cases, the ellipticals 
will be EROs at any redshift beyond ~ 0.6 in the adopted 
cosmology. Furthermore, the stellar population has to be 
formed before z ~ 1 and to be consistent with cosmology, 
an ERO-classified elliptical has to be younger than 7-8 
Gyr. 

2. The starburst galaxies can be detected as EROs at an extinc- 
tion of £(7? - V) > 1 at the very highest redshifts. When 
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E(B -V)=2 the starbursts have R - K > 5 at redshifts 
above -0.1. 

3. The 30 colour-colour combinations that can be formed 
from three of the RUHR filters exhibit the feature that 
permutations of the filters yield very similar overlap maps. 
Among the configurations formed from these filter triplets, 
the ones having the centre waveband as pivot band (i.e. the 
H band in the R - H vs. H - K configuration) have the best 
performance as discriminators. 

4. Of special interest is the R - H vs. H - K configuration. 
This provides a minimal choice of filters, as it contains the 
fundamental ERO criterion, and has a very low degree of 
overlap between the galaxy populations up to redshifts of 
2.9 for E(B - V) < 2. Using the alternative ERO criterion 
R - K > 6, the R - H vs. H - K configuration works for 

V)<2.4andz<3. 

The results in Sect. 14.21 of which item0]is a subset, repre- 
sent the limit of what is achievable with this method of photo- 
metric classifications of EROs. Since the chosen star formation 
histories are extremes (passive evolution versus constant star 
formation), the separation of the colours of the modelled galax- 
ies is maximised. A galaxy whose star formation history have 
contributions from both modes of star form ation will be more 
difficult to classify (c.f. Pierini et al. 2003). This photometric 
method is useful for galaxies clearly dominated by either star 
formation mode, but to securely establish the nature of an ERO 
other methods have to be applied as well. 
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